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The strongly coupled phase of Yang-Mills plasma with gauge group St/(3) is studied in a T- 
matrix approach. The existence of lowest-lying glueballs, interpreted as bound states of two 
transverse gluons (quasi-particles in a many-body setup), is analyzed in a non-perturbative scat- 
tering formalism with the input of lattice-QCD static potentials. The relevance of the singlet and 
the (colored) octet and 27 channels at finite temperature is discussed. We compute the equa- 
tion of state of the system in Dashen, Ma and Bernstein's formulation of statistical mechanics 
and compare to quenched SI/ (3) lattice data. Further analysis for the general case of SU(N) is 
envisaged. 
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1. Introduction 

The aim of the present work is to study the effect of two-body interactions in the thermo- 
dynamic features of the Yang-Mills (or gluon) plasma. We take recourse of Dashen, Ma and 
Bernstein's formulation of statistical mechanics ^ in terms of the 5- (or T-) matrix operator of 
gluon-gluon scattering. Correlations (leading to glueball bound states in the plasma) are studied 
within a non-perturbative approach which allows to investigate the behavior of the system in a 
range of temperatures where it is strongly interacting. We start by describing our setup for gluon- 
gluon interaction at finite temperature, and then apply our formalism to calculate the equation of 
state of the system. Next we discuss our result for the glueball spectrum in different color channels 
with special emphasis in dissociation temperatures and correlations above threshold. Finally, we 
compare our calculation of the pressure with recent lattice data in the pure gauge sector and discuss 
further applications of the present approach. 



2. Two (quasi-)gluon scattering in a T— matrix approach 

We follow Jacob and Wick's helicity formalism [Q] in order to describe a two-gluon state, the 
gluons being considered as transverse spin-1 bosons. As shown in [J3j] , two-gluon states can be 
organized in four families of helicity states in a J p basis, separated in helicity singlets |,S±;./ P ) 
and doublets \D±;J P } (symmetrization for identical two-particle states imposes selection rules on 
J). Obviously, one cannot include all the possible J p channels contributing to the dynamics of the 
system. A valid criterion to select the most relevant ones is to keep low on the averaged orbital 
angular momentum. This has been used to interpret the mass hierarchy of the glueball spectrum [^] 
(in a simple nonrelativistic picture, it provides the strength of the centrifugal barrier in scattering 
theory). With (1? ) = 2 one has, in terms of the standard \ 2S+l Lj) basis, 



S+;(F 



1/2 



\%) + 
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| 5 D () ), \S-;0-) = -\ 3 P Q ), (2.1) 



which in the singlet channel correspond to the ++ and + glueballs, respectively. With (1? 
4 one has the symmetric state |D + ;2 + ) plus three antisymmetric states: \S + ; 1 _ ), |5'_;1 + ) and 
|D_;2~). The symmetric 2 + state corresponds to the singlet 2 ++ glueball, which is also among the 
lightest ones in vacuum. We shall keep no further states in our study. 

The dynamics of gluon-gluon scattering in the plasma is encoded by the in-medium two-body 
scattering amplitude or T-matrix, which satisfies the (3D-reduced) Bethe-Salpeter equation (in 
partial waves) 

1 1'°° 

Tj P (E;q,q , )=Vj P (q,q') + j—^j o dkk 2 Vjr(q,k)G(E;k)Tjr(E;k,q'), (2.2) 

with q (q') the in-coming (out-going) momentum of a gluon in the center of mass reference frame 
and E the total energy of the gluon pair. The two-body interaction potential in momentum space, 
Vjp, is related to the potential in coordinate space, V(r), by Fourier transform and partial wave 



projection. The resolvent function in Eq. ( |2.2| ) is given by the two-gluon propagator, 



m 2 g 1 

G{E ' ,k) = e(k)e(k) 2 -E 2 /4 + 2ie(k)Z (2-3) 
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with the gluon dispersion relation e(k) = ^Jk 2 + mj and £ is related to the imaginary part (i.e., 
decay width) of the single-particle gluon selfenergy in the plasma. 



m v 4m 



(3.1) 



3. Equation of state of the Yang-Mills plasma 

The T-matrix approach is particularly suited for a straightforward implementation of the con- 
tribution of two-body interactions in the partition function of the system. At vanishing chemical 
potential, the grand potential £1 of an interacting particle gas is given by ^ 

1 

27^82 

The first term, Q. , is the grand potential of a gas of free relativistic particles (i.e., when the inter- 
actions are turned off). The second term accounts for interactions in the plasma. The sum runs 
on all the quantum numbers (v) necessary to fix a scattering channel. Below threshold, the contri- 
bution from bound states leads to the £l v term, behaving as free additional species in the plasma. 
The third term accounts for the scattering contribution above threshold. The operator enforces 
symmetrization in those channels involving identical particles, and the subscript c means that only 
connected scattering contributions are taken into account. The scattering operator is defined as 
5=1— 2ni8 (e — Ho)T, with Hq is the free Hamiltonian of the system. It is worth mentioning that 
the scattering contribution sets in at the level of the second order coefficient in a virial expansion 



of the partition function. Also note that Eq. (|3JJ) can be rewritten, by means of unitarity of the 
S-matrix, in terms of a weighted thermal average of scattering phase shifts. Once the off-shell 
T-matrix is known, the pressure is simply given by p = —Q.. Other thermodynamic observables 
can be derived form p. For example, the trace anomaly A = e — 3p and the entropy density s read, 



respectively, A= —fa [dp (P A p)] YtPll and s = -j8 2 [dpp 
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Figure 1: Left: Static free energy F\(r,T) of a quark-antiquark pair bound in a color singlet, computed 
in SU(3) quenched lattice QCD [Q], Solid lines are given by an appropriate parameterization [Q]. Right: 
r-matrix for gg scattering in the ++ singlet and octet channels. From left to right the temperatures are 
(1.05,1.10,1.15)7;, 

4. Glueballs above T r 



Following the idea in [|6|] we take recourse of lattice QCD simulations to have an estimate of the 
finite-temperature gluon-gluon interaction potential from first principles. Accurate computations 
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of the static free (Fi) and internal (U\) energies of a massive quark-antiquark pair in a color singlet 
state are available [|J], particularly in quenched St/ (3) simulations, which come handy for our 
puipose (see Fig. [j], left panel). Whether it is F\ or U\ (or none) the most appropriate potential 
to be identified with the in-medium two-body interaction at finite temperature cannot be decided 
a priori from field-theoretical arguments. Still, as in [g], we choose U\, which typically provides 
a more attractive interaction leading to larger dissociation temperatures (and thus results using 
this potential should be understood as an upper bound). We assume color scaling and derive the 
interaction potential between two quasi-gluons in color channel ^ as 

V(r,T) = ^-pi(r,T)-U 1 (oo,T)], (4.1) 

Kqq 

where K qq - is the color factor of the singlet quark-antiquark pair, and K<jf- gg = (Cf — 2N)/2N for 
SU (N), with the quadratic Casimir of the gluon pair in representation c € . As is customary, the 
potential is normalized to zero as to ensure convergence of the scattering equation. 

The only remaining unknown in our approach is the (effective) quasi-gluon mass. At finite 
temperature, it is well accepted that the infinite separation distance plateau in U\ may be interpreted 
as an in-medium contribution (selfenergy) to the gluon mass [^]. The minimal implementation is 
a shift to the bare (zero temperature) gluon mass, which at the same time has to be fixed, for 
instance, by reproducing the glueball spectrum in vacuum. In terms of the asymptotic value ofU\, 
we propose the mass shift to be 

m g (Tf = K) 2 + S(7f , 8(T) = yf^^ 1 • (4-2) 

We neglect at the moment the momentum dependence of this correction, which could be addressed 
in a full many -body study of the gluon dispersion relation. 8(T) scales with the one-body Casimir 
as suggested by the single-particle character of the correction (using two-body scaling, as in the po- 
tential, leads to a channel dependent mass correction). The structure of the color factor, ^/C g /C q , 
leads to the correct large-A^ limit of the gluon mass (as it is the case for the Hard-Thermal-Loop re- 
sult at high temperatures). Also note that this structure emerges if the gluon selfenergy is calculated 
from the T-matrix in a Brueckner scheme. 



State 


Lattice [§] 


T-matrix 


CGQCD [01 


0++ 


1.73 (5)(8) 


2.17 


1.98 


o-+ 


2.59 (4)(13) 


2.39 


2.22 


2++ 


2.40 (2.5)(12) 


2.34 


2.42 



Table 1: Masses (in GeV) of the lowest-lying glueball states at zero temperature. Our results (third column), 
are compared to the lattice data of Ref. [^] and to the Coulomb gauge QCD study in [[J. 

First of all, we test that our approach provides reasonable results in vacuum and, thus, is a 
solid stalling point for the finite temperature study. We use the "funnel" potential form with string 
breaking setting in at V s b — 2 GeV as suggested by the value of (twice) the lightest gluelump mass. 
Our results for the ground states in ++ , 0~ + and 2 ++ channels are shown in Table [j], together with 
lattice results from [[|] and the Coulomb-gauge QCD calculation in [||] (which formally shares 
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many features with our T-matrix formulation). A few comments are in order: the parameters of the 
potential, a and a, are optimized for the best fit of the lattice data for F\ beyond T c . This justifies 
the slightly too high value of the energy of the ++ glueball. A more precise agreement with the 
lattice glueball spectrum (including excited states) can be achieved if tuning the vacuum potential 
parameters independently of the finite temperature fit. More details are to be discussed in [||]. 

In Fig. p] (right panel) we show the T-matrix, at several temperatures, for the singlet and octet 
scalar glueball channels. In this plot m g is set to 1.22 GeV (namely, we use m° ~ 0.7 GeV and 8 = 
V s b/2, as in vacuum), thus the medium effects stem only from the temperature dependence of the 
potential. The goal is to illustrate the mechanism of dissociation by color screening: as temperature 
is increased, the bound states progressively move towards the two-particle threshold and eventually 
dissolve by turning to scattering states. The red, dashed curve represents det(l — VG) and vanishes 
exactly at the bound state energy. It also signals resonant states beyond threshold. 
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Figure 2: Left: ImT for gg scattering in the ++ singlet channel with in-medium quasi-gluon masses (see 
text). Right: Pressure of the SU (3) gluon-glueball gas. Lattice data from [[□]]. The dashed curve represents 
the pressure of the free quasi-gluon gas (no gg interactions). 

Next we turn the in-medium gluon mass correction on, 8(T). The bare gluon mass parameter 
is kept to m° g ~ 0.7 GeV. This is determined by our best comparison to the pure gauge lattice QCD 
data on the pressure, p = —Q. (see discussion below). The results for the imaginary part of the 
7-matrix in the singlet scalar channel are displayed in Fig. || (left panel). Two competing effects 
are responsible for the temperature evolution of the spectrum: reduction of the binding energy and 
downward shift of the threshold energy. Overall, the singlet scalar bound state experiences a mild 
shift to lower energies and dissociates at about 1.3 T c (this is the value from which detF does not 
vanish anymore; still, considerable strength remains at threshold up to about 1.5 T c ). This is in 
qualitative agreement with the spectral function analysis of Euclidean correlators by the CLQCD 



Collaboration pOfl. We also find bound states in the scalar octet channel, although considerably 
less bound. We find bound states in the pseudoscalar singlet and octet channels lying right below 
the threshold energy at the lowest considered temperature of 1.05 T c . The differences between 
singlet and octet channels are attributed to the strength of the potential, which scales with the two- 
body Casimir. In particular, note that the 27 channel is repulsive, and thus no bound states are 
found. Also note that states in the pseudoscalar channels, which in our approach correspond to 
pure P-wave scattering, are just mildly bound due to the centrifugal barrier. The tensor states, 
having an S-wave component, lie between the scalar and pseudoscalar channels regarding binding 



5 



Glueball states in the YM plasma within T -matrix approach 



Daniel Cabrera 



and dissociation temperatures. Finally, we compute the partition function of the gluon-glueball 
gas in the formalism of Dashen et al, as discussed above. The results, normalized to the Stefan- 
Boltzmann pressure, are shown in Fig. || (right panel). Note that, strictly speaking, the only free 
parameter in the calculation is the bare gluon mass, nf g . A value of 0.7 GeV provides a remarkable 



agreement with the quenched SU (3) lattice data from Panero [11]. In order to estimate the effect 
of two-body interactions, the partition function for a free gas including thermal mass corrections is 
also shown. Note that our bound states dissolve relatively quickly and only contribute significantly 
at temperatures at most up to ~ l.5T c . This illustrates the relevance of accounting for strong 
correlations from scattering in the continuum, a contribution which is readily addressed in our 
T -matrix formalism. 

5. Summary and outlook 

The relevance of gluon-gluon interactions beyond the critical temperature in the pure gauge 
SU (3) plasma has been addressed in a non-perturbative T -matrix many-body framework with the 
input of (Casimir-scaled) potentials from thermal lattice QCD. We find glueball bound states in the 
singlet and octet channels surviving up to temperatures of about 1.3 — 1.5r c , together with sizable 
threshold effects due to strong correlations beyond the two-particle threshold. With a minimal num- 
ber of parameters, we reproduce the equation of state of the gluon-glueball gas in good agreement 
with recent quenched SU (3) simulations. 

The T-matrix formalism can be applied to calculate bulk thermodynamical properties of the 
system such as the sheer viscosity, which can be easily computed in relaxation-time approximation 
within a quasi-particle picture. Our approach also allows for a straightforward extension to SU (N) 
analysis. A study of the exceptional group G2, of theoretical interest related to the origin of (de)- 
confinement, is also in progress. 
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